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Abstract 

The E x B drift motion of charged test particle dynamics in the Scrape Off Layer (SOL)is analyzed to investigate 
a transport control strategy based on Hamiltonian dynamics. We model SOL turbulence using a 2D non-linear fluid 
code based on interchange instability which was found to exhibit intermittent dynamics of the particle flux. The 
effect of a small and appropriate modification of the turbulent electric potential is studied with respect to the chaotic 
diffusion of test particle dynamics. Over a significant range in the magnitude of the turbulent electrostatic field, a 
three-fold reduction of the test particle diffusion coefficient is achieved. 
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1. Introduction 

Cross-field turbulent transport in the Scrape-Off- 
Layer (SOL) has a strong influence on divertor ef- 
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ficiency and thus on the overall performance of dis- 
charges. Recent experiments as well as numerical 
investigations suggest a non diffusive transport in 
the SOL. The transport seems to be characterized 
by large intermittent events which propagate ballis- 
tically with large velocities. On the one hand SOL 
turbulence can cross H-mode barrier and connect 
SOL and core turbulence, destabilizing the pedestal 
stable region. On the other hand, when large inter- 
mittent events propagate in the far SOL, they can 
modify the recycling pattern and consequently the 
wall particle content and thus the wall tritium inven- 
tory. Furthermore, such a convection of hot ion over- 
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densities can lead to localised heat deposition on 
components that are not designed to sustain mean- 
ingful plasma energy deposition. A strategy to con- 
trol intermittent transport can therefore prove to be 
highly valuable to achieve high performance opera- 
tion. 

Several experiments [1,2] have been dedicated to 
the control of turbulence in magnetized plasmas. Al- 
though a modelling of some of the experiments has 
been achieved, there was to date no theoretical back- 
ing for such a procedure. Our aim is to develop such 
a theoretical background and proceed by steps to- 
wards a realistic turbulence control scheme. In the 
present paper, we apply the control strategy derived 
for Hamiltonian systems to the transport of test par- 
ticles in an electrostatic potential which exhibits key 
features of the turbulent electrostatic potential gen- 
erated by the Tokam code. This is a 2D non-linear 
fluid code based on interchange instability [3,4]. 

The investigation of passive tracer dynamics has 
demonstrated to be a powerful tool for a deeper un- 
derstanding of fluid and plasma turbulence, even 
if, in the case of plasmas, one has to be cautious 
about the limits of its applicability due to the back- 
reaction of turbulence on electro-magnetic fields. 
[5,6,7,8]. 

Our strategy is to define a small apt modification 
of the system which enhances confinement, at low 
additional cost of energy [9,10,11]. In the present 
case it consists in reducing the chaotic diffusion of 
charged test particles advected by the E x B drift 
motion. This effect is obtained by adding a small 
control term to the turbulent electric potential and 
comparing the chaotic diffusion of the test particles 
advected by the original electric potential and by 
the modified one. 

In the case of turbulence with intermittent bursts 
of transport, one has to specify the control strategy: 
Either one aims at controlling turbulence in between 
the bursts of transport with a fluctuation level of 
the order of 2 %, or one aims at controlling the in- 
termittent events with fluctuation levels larger than 
10 %. The former is comparable to the case of ho- 
mogeneous turbulence and is readily achieved while 
the latter is more challenging and is presently inves- 
tigated. 

The paper is organized as follows: In Section 2 
the 2D model for SOL turbulence is recalled as well 
as its nonlinear transport properties. In Section 3 
the procedure for the computation of the control 
term is described and the results of its application 
to the E x B drift motion of charged test particles 



are presented. 

2. Interchange turbulence in the Scrape Off 
Layer 

SOL turbulence is specific with respect to core 
turbulence due to the sheath boundary conditions 
that strongly reduce the parallel current flowing on 
each flux tube. Regimes with reduced parallel cur- 
rent are then very sensitive to charge separation via 
the curvature drift since the parallel current is less 
effective to balance the charge separation. In fact, 
the interchange instability, hence curvature charge 
separation combined to a density gradient, in con- 
junction with the sheath resistivity have been pro- 
posed to explain the very large fluctuation levels 
reported in the SOL [12,13]. We consider here this 
model of SOL turbulence in the limit <C T e . In 
this model two balance equations are used, the par- 
ticle balance that governs the density transport, and 
the charge balance which takes the form of an evo- 
lution equation for the vorticity A<j>, where <f> is the 
electric potential, due to the ion inertia terms. Fur- 
ther simplifications arc introduced in the flute limit, 
hence with weak parallel gradients such that the av- 
erage along the field lines can be performed. The 
system is thus reduced to a two dimensional and two 
field model, n the density and <f) the plasma poten- 
tial normalised to T e /e. 

(J- t - DV 2 ± + {0}) n = a n n + S 

n^-fVi + {<^V3> + 00 w n = <70n (1) 

In these equations time is normalised to the ion 
Larmor frequency Qj and space is normalised to the 
so-called hybrid Larmor radius p s . D and v stand 
for the collisional transverse diffusion and viscos- 
ity (normalised to the Bohm values), S is the den- 
sity source term, here localised radially and con- 
stant in time and along the poloidal angle 9. a n and 
dtp are the sheath controlled particle flux and cur- 
rent to the wall, a n = a exp (A — <p) and u^=a{\ — 
exp (A — <p) ) where a — p s /qR stands for the nor- 
malised saturation current. The operator {<p}F = 
{(f>, F} — d x <pd y F — d y (j)d x F is the electric drift con- 
vection term. At first order in the fluctuation magni- 
tude, sheath loss terms can be linearised to yield the 
Hasegawa-Wakatani coupling term a (n — </>) [14]. 
The interested reader can refer to Refs. [3,4] for more 
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Fig. 1. Snapshot of the electric potential <f>(x, y, t) for t = 100 
obtained by Eq. (1) on a spatial grid of 512p s X 512p s . 
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Fig. 2. Typical power spectra of the electric potential <f> 
represented in Fig. 1 at t = 100 using the same parameter 
as in Fig. 1. 



details on Eqs. (1), the assumptions of this model 
and its transport properties. 

Large systems with 512p s in both the radial and 
poloidal directions are simulated for about 200 SOL 
confinement times (typically 1500 turbulence corre- 
lation times). A localised particle source, with Gaus- 
sian shape in the radial direction and width 8.5p s 
and no dependence in poloidal angle or time, splits 
the box in two regions, an unstable region where the 
average density gradient is negative and exceeds the 
threshold of the interchange turbulence, and a sta- 
ble region with positive gradient of the average den- 
sity and where the system is stable with respect to 
the interchange mechanism. A contour plot of the 
electric potential at a given time is shown in Fig. 1. 
Typical power spectra of the potential in the x and 
y direction are shown in Fig. 2. 

In the following Section we study the effect of the 
control term on E x B drift motion of charged test 
particles. 



3. Turbulent potential and test particle 
dynamics 

In the guiding center approximation, the equa- 
tions of motion of a charged test particle in presence 
of a strong toroidal magnetic field and of a time de- 
pendent electric field are [15] 

d I x\ c „ . . „ 

X =4J=^ E(X ' t)XB 

= c_f-d y 4>{x,y,t)\ 
B\d x <f>{x,y,t) J 1 [ ' 

where <j) is the electric potential, E = — and 
B = Be z . We study the 2D dynamics of test parti- 
cles obeying to Eq. (2), where the considered turbu- 
lent potential is generated by Eqs. (1). We consider 
a reduced selected spatial region and we impose pe- 
riodic boundary conditions in x, y directions. More 
precisely we choose a square box of 256,0 s x 256p s , 
not too close to the source region (in Fig. 1 the se- 
lected square box is between x £ [100, 356] in the 
radial direction and y e [1,256] in the poloidal di- 
rection). This allows one to avoid the complex anal- 
ysis of diffusive transport in finite systems, namely 
where the sources and sinks play an important role. 
With such an infinite system, one can use a straight- 
forward calculation of diffusion coefficients to char- 
acterize the turbulent transport. Since we are inter- 
ested in long time behaviour of test particle dynam- 
ics, we impose periodic boundary conditions also in 
time. Periodicity in time provides a means to gener- 
ate appropriate statistics with long time series. The 
underlying assumption is that all meaningful time 
scales of the system are shorter than the chosen pe- 
riod. A snapshot of the turbulent potential in the 
selected spatial region after the filtering procedure 
for imposing periodic boundary conditions is shown 
in Fig. 3. In what follows we denote V(x, y, t) the fil- 
tered potential. Therefore the potential V (x, y, t) is 
L x = 256p s periodic in the x-direction, L y = 256p s 
periodic in the y-direction and T — 25650 fi" 1 pe- 
riodic in time. 

In order to have a 1-pcriodic dependence in 
(x,y,t) for simplicity, we perform the change of 
variables 

x y ~ t , , 

x = T, v = j-i t= T- (3) 

±Jy 1 

It follows that 

dx _ T dx _ T dV 

~dl ~ T x ~di ~ ~T x l)y~' 
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Fig. 3. Snapshot of the electric potential V(x, y, t) for t = 100 
over the spatial grid of 256p s X 256p s obtained from the 
original potential </> represented in Fig. 1 after the filtering 
in space and time. 



From 

dV _ dVdy 
dy dy dy' 

we get 



dx 
dt 



T dV(x,y,t) 



dy 



(5) 



(6) 



Thus the equations in the (x,y, t) variables are 
equivalent to the equations in the (x, y, t) variables 
if the potential is multiplied by a factor L T L . 

The test particle dynamics governed by Eq. (2) is 
studied with the aid of numerical simulations. Since 
the potential is given over a spatio-temporal grid, 
the estimation of the electric field which acts on the 
particle trajectories is computed using a four point 
bivariate interpolation in space and a linear inter- 
polation in time [16]. The time integration is per- 
formed using a fourth order Runge-Kutta scheme. 

We consider a set of A4 particles (A4 of order 
1000) uniformly distributed at random at t = 
and we study the diffusion properties of the system. 
We verified that for the values of the parameters in 
the range of interest the mean square displacement 
{r 2 {t)) is a linear function of time and we computed 
the corresponding diffusion coefficient 



D* 



lim ■ 

t— >oo 



2 W> 



t 



as a function of the amplitude of the turbulent po- 
tential. In Fig. 4 the results obtained for small values 
of the electric potential are shown. 

The first step in the computation of the appropri- 
ate modification of the potential which entails a re- 



Fig. 4. Test particle diffusion coefficient D* versus a (mag- 
nitude of the potential) in log-lin scale obtained for the tur- 
bulent potential V without (full squares, solid line) and with 
(full circles, dashed line) the addition of the control term J2 
given by Eq. (10). 

duction of chaotic diffusion is the Fourier expansion 
of the electric potential, i.e. 

V(x,y,t)= V klk2k3 e 2 ^ kix+k ^ +k3t \ (7) 

fci ,k2 ,k3 

Then we compute fV, that is the time- integral of V, 



TV 



E 



Vk 1 k 2 k 3 2wi(kix+k 2 y+k 3 t) 



(8) 



and the constant part of V with respect to time that 
is written 1ZV: 



fei ,fe 2 



2-rri(k-ix+k2v) 



(9) 



The control term / 2 is given by the Poisson bracket 
between TV and (K + 1)V, 



--{TV,(U+1)V} 

1 fdTV d{Tl + l)V dTV d(K + T)V 



2 V dx 



dy 



dy 



dx 



(10) 



A contour plot of at a fixed time is depicted in 
Fig. 5. Finally we study the test particle dynam- 
ics in the modified potential V + fa. The compari- 
son between the diffusion coefficient obtained with 
and without the control term is shown in Fig. 4. 
One readily observes two main features, on the one 
hand the control term exhibits higher frequencies 
than the turbulent electrostatic field, typically a fac- 
tor 2 larger has expected from the quadratic form 
of the control term, Eq.(lO), on the other hand the 
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Fig. 5. Contour plot of the control term f2(x, y, i) for t = 100 
given by Eq. (10). 

magnitude of the control term is more than a factor 
10 smaller. This reduction compensates for the in- 
crease of the electric field via the increased frequen- 
cies. A significant reduction of D* is observed in the 
controlled case. These preliminary results obtained 
for small amplitudes of the electric potential show 
that the control strategy we propose can be applied 
to much more turbulent flows than the one consid- 
ered in [9]. Further investigations on the efficiency 
domain of the control term and its robustness are 
required. 

4. Discussion and conclusion 



References 

[I] C. Schroder, T. Klinger, D. Block, A. Piel, G. 
Bonhomme, V. Naulin, Phys. Rev. Lett. 86, 5711 (2001) 

[2] J. Stockel, et al, Plasma Phys. Control. Fusion 47, 635 
(2005) 

[3] Y. Sarazin and Ph. Ghendrih, Phys. Plasmas 5, 4214 
(1998) 

[4] Ph. Ghendrih et al, Nucl. Fusion 43, 1013 (2003) 

[5] V. Naulin, O.E. Garcia, M. Pricgo, J. Juul Rasmusscn, 
Phys. Scr.T122, 129 (2006). 

[6] V. Naulin, Phys. Rev. E 71, 015402 (2005). 

[7] M. Priego, O.E. Garcia, V. Naulin, J. Juul Rasmusscn, 
Phys. Plasmas 12, 062312 (2005). 

[8] R. Basu, V. Naulin, J. Juul Rasmussen, Commun. 
Nonlinear Sci. Num. Simul. 8, 477 (2003). 

[9] G. Ciraolo et al, Phys. Rev. E 69, 056213 (2004). 

[10] C. Chandre et al, Phys. Rev. Lett. 94, 074101 (2005). 

[II] C. Chandre, M. Vittot, G. Ciraolo, Ph. Ghendrih, R. 
Lima, Nuclear Fusion, 46, 33 (2006) 

[12] A.V. Ncdospasov, Sov. J. Plasma Phys. 15, 659 (2003). 

[13] X. Garbct et al, Nucl. Fusion 31, 967 (1991). 

[14] M. Wakatani and A. Hasegawa, Phys. Fluids 27, 611 
(1984). 

[15] T. J. Northrop, Annals of Physics 15, 79 (1961). 

[16] M. Abramowitz and LA. Stegun, Handbook of 
Mathematical functions, Dover Publications, (1970). 



In this paper, we have shown that the procedure 
developed to control turbulent diffusive transport 
of test particles can be applied to particles trans- 
ported in a realistic turbulent electrostatic poten- 
tial. The procedure to compute a low magnitude 
control term provides an efficient turbulence reduc- 
tion is presented. For this test of the control strat- 
egy the turbulent field is simplified to simplify the 
analysis of test particle transport. More generally 
the proof of principle that a control strategy can be 
efficient is tested on the simplest available model. 
When comparing the diffusion coefficient of the test 
particles in the uncontrolled system to that in the 
controlled typical factor 3 reduction of the 

diffusion coefficient is achieved with a factor 10 lower 
magnitude of the control term. Far more investiga- 
tion is still required on the way towards realistic con- 
trol systems for the turbulent transport, however, 
the present paper shows that a theory based start- 
ing point is available. 
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